GLOBAL SOLUTIONS TO BUBBLE GROWTH IN 
POROUS MEDIA 



LAVI KARP 

Abstract. We study a moving boundary problem modeling an 
injected fluid into another viscous fluid. The viscous fluid is with- 
drawn at infinity and governed by Darcy's law. We present solu- 
tions to the free boundary problem in terms of time-derivative of a 
generalized Newtonian potentials of the characteristic function of 
the bubble. This enables us to show that the bubble occupies the 
entire space as the time tends to infinity if and only if the inter- 
nal generalized Newtonian potential of the initial bubble is a qua- 
dratic polynomial. Howison [7], and DiBenedetto and Friedman 
[2], studied such behavior, but for bounded bubbles. We extend 
their results to unbounded bubbles. 



1. Introduction 

The propose of this note is to apply the technique of generalized 
Newtonian potential to the study of bubble growth in porous media 
and Hele-Shaw flows. This technique enables the computation a po- 
tential of measures with unbounded support in a similar manner to 
the ordinary Newtonian potential . It is a multi-valued right inverse 
to the Laplacain and unique up to a harmonic polynomial of degree 
not exceeding two (see [8, 10, 13]). We can therefore compute the po- 
tential of the characteristic function of arbitrary large set IR n . The 
Newtonian potential theory is a basic tool in the studies of Hele-Shaw 
flows [5, 6, 14, 18, 21], and since in a bubble growth in porous media 
the fluid motion is in unbounded domains, the application of the gen- 
eralized Newtonian potential is rather natural in this type of moving 
boundary problems. 

Consider ~R n as homogeneous porous medium filled with a viscous 
fluid. Another fluid is injected and forms a bubble which occupies a 
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domain D(t) at each time t. The fluid withdrawn at infinity at a certain 
rate and the bubble D(t) increases with the time t. Let Q(t) = M. n \D(t) 
and \&(x, t) denotes the pressure of the incompressible fluid in Q(t). We 
will prove that ^>(x,t), the solution to this moving boundary problem 
(1.1) below, is a time-derivative of a generalized Newtonian potential 
of XD(t)i the characteristic function of bubble D(t). In this formulation 
it does not matter whether the bubble is bounded or unbounded. 

The representation of solutions by potentials has several advantages 
specially in higher dimensions where the tool of conformal mappings 
is not available. For example, this enables us to construct in a simple 
manner solutions for which the bubble exists for all for all t > and 
occupies the entire space as as t — > oo. This problem was settled by 
Howison [7], and DiBenedetto and Friedman [2], but under the condi- 
tion that the initial bubble is a bonded domain. Here we show that this 
type of fluids motions exists if and only if the internal generalized New- 
tonian potential of the initial bubble Xd(o) is a quadratic polynomial. 
It thus extends their result to unbounded bubbles. 

The known examples of domains in M n for which internal generalized 
Newtonian potential of their characteristic set equals to a quadratic 
polynomial are (a) ellipsoids, (b) convex domains bounded by ellip- 
tic paraboloids, (c) domains bounded by two parallel hyperplanes, (d) 
cylinders over (a) and (b), and (e) half-spaces [10]. The complements 
of these domains are null quadrature domains, that is, an open set O 
for which 

hdx = 

for all harmonic and integrable functions h in Q. In the two dimen- 
sional plane Sakai classified those domains [17], but the classification 
in higher dimensional spaces is an open problem. Dive [4], and Nikli- 
borc [15], proved that if the internal Newtonian potential of a bounded 
domain is a quadratic polynomial, then it must be an ellipsoid (their 
proof in given in M 3 , for proofs in arbitrary dimension see [2, 9]). The 
classification problem of this type of domains is also settled under the 
a priori assumption that the domain is contained in a cylinder of co- 
dimension two [10]. A recent progress in the classification problem was 
obtained in [11] and for its current state see Remark 4.3 below. 

1.1. The formulation of the moving boundary problem. The 

bubble is formed by injected a fluid of negligible viscosity into another 
incompressible viscous fluid. We denote the set occupied by the bub- 
ble at time t by D{t) and fi(i) = R n \ D(t). Suppose W 1 to consists 
of homogeneous porous medium, then the motion of the flow in Q(t) is 
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subject to the Darcy's law which asserts that the velocity v is propor- 
tional to the gradient of the pressure 



v = -kVj*. 
Since the flow is incompressible, div v = 0, hence 

= in fi(i) for all t > 0. 

Assuming there is no surface tension, then the pressure is constant on 
the boundary dQ(t), so we may set 

V(x, t) = x G D(t) for all t > 0. 

The free boundary moves with velocity — V x \l/, that is 

— = -V„ on 

where rj is the outward normal pointing into Q(t) and is is the 
outward velocity of the free boundary dD{t). If D(t) = {g(x,t) < 0}, 

then f£ = • f -r^^-^) and V„ = ^H. Generalized Newtonian 

potential of densities in L°°(M ri ) have |x| 2 log|x| growth as x — > oo 
[8, 20], and since in our setting the bubble D(t) may has infinite volume 
we allow ty(x,t) = o(|x| 3 ) at infinity. We summarize these conditions: 



(1.1a) A# = in Q(t), 

(1.1b) ^ = in D(t), 

(i.ic) — = -v„ on acitt), 

Ofj 

(l.ld) *(x, t) = o(|x| 3 ) as \x\ -> oo, 

(l.le) V^(x,t) = o(|x| 2 ) as |x| ->• oo. 



2. Generalized Newtonian Potential 

The definition and establishment of basic properties of this potential 
were carried out in [10]. Here we recall the definition and present an 
estimate which is needed for the current application. 

The Newtonian potential of a measure fi with compact support is 
defined by means of the convolution 



(2.1) V(jj)(x) = (J * fj) (x) = / J(x - y)dii(y) 
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where 

{-— \og\x\ , n = 2 

27T 1 
7 s r~, — n i n > 3 
(n - 2)w n |x| ri - 2 

and oj n is the area of the unit sphere in IR n . The potential V(fi) sat- 
isfies the Poisson equation AV(fi) = —\x in the distributional sense. 
Similarly, the generalized Newtonian potential is a right inverse of the 
Laplacian, but it is multi-valued and acting on the space £, the space 
of all Radon measures /i in M. n satisfying condition 

( 2 - 2 ) Mc-=J YTJiW 

Note that C contains all measures with densities in L°°(IR n ). The linear 
space C is the Banach space with the norm defined by (2.2). For 
fi G C we first define V a (n), the operator of the third order generalized 
derivatives of the potential: 

(2.3) <p):=-J d a V(<p)(x)dfx(x), tpeS, \a\ = 3. 

Here S is the Schwartz class of rapidly decreasing functions and V a (fi) 
is a tempered distribution. 

Definition 2.1. The generalized Newtonian potential V\jj] of a measure 
/i G £ is the set of all solutions to the system 

Am = —fi 

d a u = V a (fi) , \a\ = 3 . 

The existence of solutions to system (2.4) was proved in [10] and it 
is unique modulo H2, the space of all harmonic polynomials of degree 
at most two. The operator V : C — > S' /H2 is continuous [10]. 

We will not distinguish here between ordinary and generalized po- 
tentials. By a potential of a set A we mean the potential of its charac- 
teristic function \a- 

Let k be a positive integer, then 

ll^llfe = SUp SUp ((1 + \x\f \d k (f(x)\) 
{\a\<k} M" V J 

is a semi-norm on S. We shall need the following estimate. 

Proposition 2.2. Let a and [3 be multi-indexes such that \a\ = 3. Then 
for any ip G S 

(2.5) \&*^v{ v ){x) < c (1 + \x\y {n+1+m yh (n+1+m , 

where the constant C does not depend on if. 



(2.4) 
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For (3 = this was proved in [10, §1]. Since only a slightly mod- 
ification is needed in order to extend it for (2.5) we leave it to the 
reader. 



3. The Generalized Newtonian Potential of the Bubble 

In this section we show that solutions to the moving boundary prob- 
lem (1.1) can be expressed in terms of the generalized potential. 

Proposition 3.1. Let V a (xD(t)) be the third order derivatives operator 
defined in (2.3) and assume \l/ satisfies (1.1), then 

(3.1) Jt Va{xD ^ = ^ H = 3 ' 

where the identity (3.1) is in the distributional sense. 

Proof. Let ip e S, then Proposition 2.2 implies that dV a (tp) E L 1 (M. n ) 
and therefore V a (xR»)(<p) = ~ f R » &*V{tp)dx = 0. Since fi(t) = E n \ 
D(t), 

(3.2) (V a ( X D(t)),<p) = [ d a V(ip)(x)dx = - [ d a V(p)(x)dx 

J D(t) Jn(t) 

and 

(3.3) 1 {(V a ( X D(t))^)) = - [ d«V{y){x)V v dS. 

Let Br be a ball of radius R and apply the Green's formula, then 



A (d a V(ip)) Vdx 

Q(t)nBjj 

By the estimate (2.5) and the growth assumptions (l.ld) and (Lie) of 
the second term of the right hand site of (3.4) tends to zero as R goes 
to infinity. Thus, letting R — > oo and using the boundary conditions 
of (1.1) we get 

(3.5) f A(d a V(<p))^dx= [ d a V{p)Vr,dS. 

Jn(t) JdQ(t) 
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Since AV(<p) = —tp for any test function, we see from (3.3) and (3.5) 
that 

(3.6) i. ((V a (xD {t) ), f)) = - [ A (d a V(ip)) Vdx = - [ dPtpVdx 
at Jn(t) J 

which is the distributional meaning of (3.3). □ 

Corollary 3.2. If the bubble D(t) occupies the entire space as t tends 
to infinity, that is, 

(3.7) lim D(t) = R n , 

t— >oo 

then the potential of the initial bubble D(0) coincides with a quadratic 
polynomial on D(0). 

Proof. Taking the integral of (3.1) we get 

(3.8) (V a ( X D(t))M - (V a (XD { o)), = ] &*<P(?Mx, s)dxds. 

Since ^f(x, s) — for all x G D(0) and all s > 0, we get from (3.8) that 

(3.9) (V a ( X D(o)), ip) = (V a (XD(t)), <p) for all <p G C?(D(0)). 

Now if (3.7) holds, then by the continuity of the generalized potentials 
we have that 

(V a ( X D { o)), f) = lim (V a ( Xm ), f) = (V a ( X ^), <p) = 

for all ip G C^°(D(0)). Hence the third order derivatives for any u G 
V[Xu(o)] vanish in D(0). □ 

The following formula is the analogous of Richardson's theorem [16], 
but for suction at infinity. This was previously proved by Entov and 
Etingof [5] (see also [21, §4.6]), in the plane and for a bounded bubble. 
We thus generalized their result to the space M n and for unbounded 
bubbles. 

Theorem 3.3. Suppose the solution of (1.1) exists fort G [0, T]), then 
there exists a generalized Newtonian potential u(-,t) G V^[xD(t)] such 
that 

(3.10) ^u(x,t) = a(t) + ^(x,t). 

Proof. Let w G V[xo(t)] an d |«| = 3, then by Proposition 3.1 

(3.11) j t d a w(x,t) =d a ^(x,t), 
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which implies 

(3.12) d a w(x,t) = F a (x) + f d a ty{x,s)ds. 

Jo 

Since diF a = djFp for any a, (5 with did a = djd@, there is a function F 
such that d a F = F a . Hence 

(3.13) w(x,t) = F(x) + I ty(x 1 s)ds + q(x 1 t), 

Jo 

where q(x,t) is a quadratic polynomial. We can write it in the form 

q(x, t) = \x\ 2 M{t) + h(x, t) + A(t), 

where h(x, t) is harmonic polynomial of degree not exceeding two (see 
e.g. [1]). Now ^(x,t) = for x E D(0) and t E [0,T], therefore in 
-D(O) we have 

(3.14) - 1 = Aw(x, t) = AF(x) + 2nM(t). 
Differentiating with respect to t gives 

(3.15) = ^-Aw(x, t) = 2n^-M(t). 

So M(t) does not depend on t. Letting u(x, t) = w(x, t) — h(x, t), then 
u belongs to V^[xd(o] anc l satisfies (3.10) with a(t) = 4:A(t). □ 

We shall now see that the converse statement to Theorem 3.3 is also 
true. 

Theorem 3.4. Let D(t) be a continuous family of domains such that 
D(ti) C D(t 2 ) for t\ < t 2 and dD(t) is sufficiently smooth so that the 
Green's identity holds. If for t E [0,T] there are potentials u(-,t) E 
V[xD(t)} such that 

(3.16) m(x, t 2 ) — u(x, ti) does not depend on x for x E D(ti), 
then the solution to (1-1) is given by 

(3.17) ^(x, t) := ^-u(x, t) - ait). 

dt 

Proof. Here we argue similar to [2]. Let (p E C^°(M. n ) and t\ < t 2 . Then 
lim / (pdx = I V v dS, 

«2->ti t 2 — ti JD{t 2 )\D(tt) JdDih) 

where V v is is the outward velocity of the free boundary dD{ti). Since 
- 1 - [ <pdx=- — l — [ (u(x,t 2 ) -u(a;,ti)) Atpdx, 

t2 — ti JD{t 2 )\D{t x ) t 2 -txJ 
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4:u(x,t) exists and 

(3.18) (^u(;t),A<p) = -{V n7 (p). 

By condition (3.16), u(x, t)—u(x, 0) does not depend on x for x G -D(O). 
So we denote u(x, t) — u(x, 0) by A(t) and define 

%{x,t) :=ju{x,t)-A!(t). 

We shall now verify that this $(x, t) satisfies (1.1). Since u(x,t) G 
V[xi3(t)], h is harmonic in Q(t) and therefore (1.1a) holds. Condition 
(3.16) implies that ^f(x, t) = for x G D(t), so (1.1b) holds too and 
the growth properties (l.ld) and (Lie) are consequence of the known 
estimates of the generalized potential of L°°(K n )-densities [8, 20]. In 
order to verify (1.1c) we use (1.1b), (3.18) and the Green's formula, 
these yield 

d f d 

-(Vr,,cp) = (—«(•, t), Aip) = / —u(x,t)A(f(x)dx 



(3.19) 



>dD{t) 

This completes the proof. □ 



4. Construction of global solutions 

In this section we use Theorem 3.4 in order to construct solutions to 
the moving boundary problem (1.1) such that the bubble D(t) will oc- 
cupy the entire space as t goes to infinity. By Corollary 3.2 a necessary 
condition for that is 

u(x) = q(x) for xE D(0), uEV[xd(o)], 

where q is a quadratic polynomial. Here we shall prove that this is also 
a sufficient condition. 

We first recall few known facts. A classical theorem of Newton says 
that the gravitational force of a homogeneous shell between two similar 
ellipsoids vanishes in the cavity of the shell (see e.g. [3, 12]). This is 
equivalent to the property that the ordinary internal Newtonian po- 
tential of ellipsoids is a quadratic polynomial (see e.g. [2]). We also 
know that the following unbounded domains have their internal gener- 
alized potential coincides with a quadratic polynomial: (i) strips, (m) 
half-spaces, (iii) convex domains bounded by elliptical paraboloid, and 
(iv) cylinders over ellipsoids and over these domains [10, 19]. 
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The first example shows that there are many motions of bubbles 
which fill the entire space as t runs to infinity. 

Example 4.1. Let E(t) be any continuous family of ellipsoids such that 
E(ti) C E(t 2 ) whenever t x < t 2 and lim*-** E(t) = R n . Then the 
Newtonian potential V(xE(t))( x ) = q(x,i) for x G E(t), where q is a 
quadratic polynomial satisfying Aq = — 1. So we may write 

n 1 

(4.1) q{x,t)=^2a i {t)x 2 i +p 2 {x,t), ai {t) + ■ ■ ■ + a n {t) = --, 

i=l 

where p 2 is a harmonic polynomial of degree not exceeding two. Then 

n-l 

h(x, t) :- 



^Oi(t)z? + (a n {t) + -j x\ 



is harmonic polynomial of degree two and therefore 

u(x, t) = V(xE(t))(x) - h(x, t) - p 2 (x, t) 

is a generalized Newton potential of XE(t) satisfying u(x, t) = —\x 2 n on 
E(t). Therefore, condition (3.16) holds so we conclude that ^>(x,t) = 
■^u(x,t) is a global solution to the moving boundary problem (1.1). 

Of course these solutions may not have a meaningful physical inter- 
pretation, since they have a quadratic growth at infinity. It is therefore 
reasonable to look for solutions which have a minimal growth at infin- 
ity. We see from (4.1) that a necessary condition for that is that the 
coefficients cij(£) are independent of t. Such behavior of the coefficients 
a,i(t) occurs when E(t) are concentric ellipsoids having there center at 
the origin (see e.g. [15]). 

Thus the ellipsoids of the form 

= eR " : fe) 2 + "' + fe) 2<(1 + (): 

will form the bubble for which ^(x, t) will have a minimal growth at in- 
finity. This example where also obtained by Friedman and DiBenedetto 
[2]. In addition, any cylinder 

\otiJ \a k J 

will also gives global solutions which has minimal a growth in directions 
off the cylinder. 

We now turn showing that elliptical paraboloid also form a global 
solution. 
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Example 4.2. Let 

(4.2) F=tx: x n > 

then any generalized Newton potential of xf coincides with a quadratic 
polynomial on F. We could therefore argue like as in Example 4.1, but 
that will not have a minimal growth. We claim that there is u G V[xi?] 
such that u(x) does not depend on x n for x G F. Assuming it for the 
moment, and define the bubble F(t) by a translation in the Xyi axis, 
that is, 

F(t) = {x G R n : {x x , z n _i, x n + t)ED}. 

Then obviously u(x, t) := u(a;i, ...,x n -x,x n +t) is a generalized potential 
of XF(t) an d since w(x) does not depend on x n , u(x,t) = u(x) for 
x G F(t). Thus u(x,t) satisfies condition (3.16) and by Theorem 3.4 
■^u(x,t) is a global solution to (1.1). In addition, lim^oo F(t) = M. n . 

In order to prove the claim, we adopt the following argument from 
[11, §4]. There we showed that if V"[xd] coincides with a quadratic 
polynomial in D and 

4-3 Km V ,, R / = 0, 

then there is a potential w G V[%£)], such that 

n n 

w(x) = q(x) := ^ biXi + c, for x G D 

ij=l i=l 

and 



(4.4) lim W{PX) ~ q{px) = f 



CtijXiXj . 



Here Vol(-) denotes the Lebesgue measure in M. n and B p a ball with 
radius p. 

Now the set F given by (4.2) certainly satisfies condition (4.3). So 
let w G V[xf] be the potential which satisfies (4.4). Since w(x) — q(x) 
vanishes on F, w(px) — q(px) vanishes on the positive for every 

p and therefore the limit in (4.4) also vanishes there. Hence a n = 0. 
Letting u{x) = w(x) — Yll=i{ a in + CL n i)xiX n — b n x n gives the required 
potential which also provides a solution with minimal growth at the 
direction of the negative 

Remark 4.3. The complete characterization of domains D with internal 
potential equal to a quadratic polynomial is not known. We recently 
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proved in [11] that if D is such domain, then it must be one of the 
following type: 



!. If iimm „ V»(f nB„ > o, then D is half-space; 



2. If Hindoo — — — — = 0, then there are two possibilities: 

2a. If D is contained between two parallel hyperplanes, then 
D is an ellipsoid, or a cylinder over an ellipsoid; 

2b. If D is not contained between two parallel hyperplanes, 
then D = {x n > f(xi,...,Xk), 1 < k < n} and / is real 
analytic convex function. 

It is obvious that we can extend Example 4.2 to domains of type 2b. 
above. We thus proved: 

Theorem 4.4. Let ty(x,t) be a solution to the moving boundary prob- 
lem (1.1). Then the bubble D(t) occupies the entire space as t goes to 
infinity if and only if the generalized Newtonian potential of the initial 
bubble V[xd(o)] coincides with a quadratic polynomial in D(0). 

Finally, these techniques can be used also for contractions of flows. 
For example, we can construct solution to (1.1) so that the paraboloid 
D in (4.2) will contract to a half-line. Here we use again the fact that 
there is u G V^[xd] that does not depend on x n for x e D. 
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